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Abstract 

We consider transformation maps on the space of states which are symmetries 
in the sense of Wigner. By virtue of the convex nature of the space of states, 
the set of these maps has a convex structure. We investigate the possibility of 
a complete characterization of extreme maps of this convex body to be able to 
contribute to the classification of positive maps. Our study provides a variant of 
Wigner's theorem originally proved for ray transformations in Hilbert spaces. 
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1 Introduction 

Symmetries play a very important role in physics, as it has been stressed by Wigner 
on several occasions [H [2j [3] . The way symmetries are realized depends on the theory 
under consideration and more specifically, according to Felix Klein, on the corresponding 
geometric structure of the carrier space, these are 'kinematical' rather than 'dynamical' 
symmetries. It is well known that any description of physical systems requires the 
consideration of states and observables along with a pairing among them providing a 
real number with a computable probability |4j. In the Schrodinger-Dirac description 
of Quantum Mechanics one associates a Hilbert space with a quantum system, states 
are identified with rays of this space and observables are a derived concept - they 
are identified with self-adjoint operators, while symmetries are defined to be bijections 
among rays which preserve probability transitions. 

In the C*-algebraic approach to Quantum Mechanics, originated from the Heisenberg 
picture, observables are identified with real elements of this C*-algebra, while states 
are a derived concept - they are identified with positive normalized functionals on the 
space of observables. The space of observables carries the structure of a Jordan alge- 
bra and this was the point of view of Kadison to define symmetries as Jordan-algebra 
isomorphisms [51 Ej. C* -algebras are quite convenient to deal with the description of 
composite systems, the dual space of states turns out to have a rather involved geomet- 
rical structure. In particular, to take into account the distinction between separable 
and entangled states on the space of states, one is obliged to give up linear superpo- 
sition in favor of convex combinations. This change of perspective introduces highly 
nontrivial problems, specific to the 'convex setting'. As shown elsewhere [TJ, in finite 
dimensions the space of states turns out to be a stratified manifold with faces of various 
dimensions. 

The aim of this paper is to deal with symmetries as those transformations on the space of 
states which are appropriate for its geometrical structure. In doing this we end up with 
yet another variant of the celebrated Wigner's theorem on the realization of symmetries 
as unitary or antiunitary transformations on the Hilbert space. The literature on this 
theorem, which is also available on text books [HI [H] in addition to the famous book by 
Wigner [10J, is huge. We limit ourselves to a partial list trying to give a sampling of the 
various approaches which have been taken in the years [TH [TJl [13], [15j [T6l [T7J, HH [19] . 

The paper is organized in the following way: in Section [2] we give a short geometrical 
description of the set of density states in a finite-dimensional Hilbert space. Density 
states form a convex body in the space of Hermitian operators. The set of affine maps 
which map a convex set K into itself, called simply positive maps, is also a convex 
set in the space of affine maps. Characterization of positive maps, e.g. by identifying 
the extremal ones, i.e. maps which can not be decomposed into a nontrivial convex 
combination of other positive maps, can lead to a useful description of the underlying 
set K. Finding extreme points of such maps is, however, a difficult task, even if we 
know explicitly the extreme points of K. In Section [3] we discuss and give examples 
of positive maps for which their extremality can be established upon analyzing the 
number of extreme points in the image. In Section H] we connect the obtained results to 
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the Wigner's theorem expressed in terms of positive maps which are bijective on pure 
states. Sections [5] and [6] are devoted to completely positive maps. In particular, we 
show again how the number of extreme points in their image establishes their form and 
extremality. We conclude with Section [7] containing illustrative examples of extreme 
positive and completely positive maps in low dimensions. 



2 Density states 

Let TC be a finite-dimensional Hilbert space, dim7i = n, and let gl{H) be the space of 
complex linear operators on 7i. The space gl(H) is canonically a Hilbert space itself 
with the Hermitian product (A, B) = Tr(A' o B). As in [7j, we shall treat the real 
linear space of Hermitian operators on 7i as the dual space, u*(7i), of the Lie algebra 
(of antihermitian operators), u(H), of the unitary group U(H). We have the obvious 
decomposition gl(H) = u(H) Q)u* (H) into real subspaces with a natural pairing between 
u(H) and u*(H) given by 

(1) (A, T)* = i- Tr(AT), (A, T) G u*(H) x u(H), 

and a scalar product induced on u*(7i) by the Hermitian product and given by 

(2) (A,B)*=Tr(AB), A,Beu*(H). 

We denote with || • ||* the corresponding norm. 
The coadjoint action of U(TL) on u*(TC) reads 

(3) A\->UAlfl, Aeu*(H), UeU(H). 



We denote by V(TC) the space of positive semi-definite operators from gl(7i), i.e. of 
those p G gl{Bi) which can be written in the form p = T^T for a certain T G gl(TC)- It 
is a cone, since it is invariant with respect to the homotheties by A with A > 0. The set 
of density states V(TL) is distinguished in the cone V(TC) by the equation Tr(p) = 1, so 
we will regard V(H) and V{7i) as embedded in u*(H). 

The space T>{H) is a convex set in the affine hyperplane u*(7i) in u*(7i), determined 
by the equation Tr(A) = 1. The model vector space for u\{H) is therefore canonically 
identified with the space of Hermitian operators with trace 0. The space 2t(u^(7i)) of 
affine maps of u\(1-i) can be canonically identified with the space of these linear maps 
$ G C(u*(7i)) which preserve the trace. 

It is known that the set of extreme points of T>{7i) coincides with the set T> X ^H) of 
pure states, i.e. the set of one- dimensional orthogonal projectors and that every 

element of T>{7i) is a convex combination of points from T> X ^H). The space T> 1 {7{) of 
pure states can be identified with the complex projective space P7i ~ CP n_1 via the 
projection 

H\{0}3x^ l -^f}-eV 1 (n), 
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which identifies the points of the orbits of the C \ {0}-group action by complex homo- 
theties. 

If we choose an orthonormal basis ex, . . . , e n in TC, we can identify u*(7i) with the real 
vector space u* (n) of Hermitian nxn matrices, u\ (7i) with the affine space of Hermitian 
n x n matrices with trace 1, U(7i) with the group U(n) of unitary matrices, T>{7i) with 
T>(n) - the convex body of density nxn matrices, etc. Recall that the dimension of 
u\{n) is n 2 — 1 and the dimension of u*(n) is n 2 . 

Almost all above can be repeated in the case when 7i is infinite-dimensional if we 
assume that all the operators in question, i.e. operators from gl(H) and u*{7i) are 
Hilbert- Schmidt operators (see |2Q]). The positive semi-definite operators then, being 
of the form AA\ are trace-class (nuclear) operators, so density states are trace-class 
operators with trace 1. There are some obvious minor differences with respect to finite 
dimensions: for instance, the convex set T>{7i) of density states is the closed convex 
hull of the set V 1 (TC) of pure states, rather than just the convex hull, etc. 

3 Positive maps of convex sets 

If K is a convex set in a locally convex topological vector space E, then the set Pos(fT) 
of those continuous linear maps $ : E — *> E which map K into K is a convex set in the 
(real) vector space C(E) of all continuous linear maps from E into E. We will refer to 
elements of Pos(K) as to linear K-positive maps, or simply to linear positive maps, if 
K is determined. 

If K is compact, then, due to the Krein-Milman Theorem, it is the closed convex hull 
of the set K° of its extreme points (points which are not interior points of intervals 
included in K), K = con(K°). In this sense, compact convex sets K are completely 
determined by their extreme points. 

However, it should be made clear from the beginning that the concepts of convex set, 
positive map, etc., are taken from the affine rather than linear algebra and geometry. 
In an affine space E, one can subtract points, x = p — p', to get vectors of the model 
vector space E = v(E), or add a vector to a point, p = p' + x, to get another point, but 
there is no distinguished point that serves as the origin. More generally, in affine spaces 
we can take affine combinations of points, i.e. combinations Yli^iPi sucn that J^Aj = 
1. If all Aj are non- negative the corresponding affine combination is just a convex 
combination. We say that points p , . . . ,p r G E are affinely independent if none is an 
affine combination of the others. This is the same as to say that p\ — po, . . . , p r — po G E 
are linearly independent vectors. 

Convex sets in our approach will live in affine spaces. In this sense the Krein-Milman 
Theorem tells us something about compact convex sets in affine spaces modeled on 
locally convex linear spaces. 

One can think that the problem is artificial, since by choosing a point in an affine space 
as the origin we end up in the model vector space. However, choosing a point is an 
additional information put into the scheme which changes our setting. The situation is 
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like in a gauge theory, where we can fix a gauge. But a fixed gauge has, in general, no 
physical interpretation, so we rather try to use gauge-invariant objects. 

The second instance of affine space presence is the fact that in many situations, even 
when we work in a true linear space, it makes much more sense to admit that positive 
maps are affine. Note that affine maps on an affine hyper space E of a linear space E 
come exactly from linear maps in E which preserve E. On the other hand, every affine 
space (or even affine bundle) E can be canonically embedded in a linear space (vector 
bundle) E as an affine hyperspace (affine hyperbundle). We refer to [21, [22] for the 
corresponding theory with interesting applications to frame-independent formulations 
of some problems in Analytical Mechanics. 

Definition 1. (a) Let Ej be a real affine space modeled on a locally convex topological 
real vector space Ei = v(Ej), i — 1,2. We say that a map $ : Ei — ► E2 is an affine 
map if there is a continuous linear map v ($) : E\ — > Ei such that for any p G Ei and 
any x G E\, we have $(p + x) — $(p) + v (<&)(x), where p t— > p + x is the natural action 
of Ei on Ei. The space of affine maps from Ei to E2 will be denoted by 2l(Ei, E2). If 
Ei = E 2 = E, for the space of affine maps on E, i.e. for 2l(E, E) we will write shortly 
21(E). 

(b) Let 21(E) be the space of all affine maps on E and let K be a convex set in E. 
By positive maps on K (or simply positive maps if there is no ambiguity about K) we 
understand these affine maps $ G 21(E) which map K into K. The set of all positive 
maps on K will be denoted by ?$(K). 

(c) By a convex body we will understand a compact convex set K with non-empty 
interior in a finite-dimensional Euclidean affine space E. 

Note that the set T>{H) of density states for quantum systems with a finite number of 
levels is an example of a convex body, as it is canonically embedded in the Euclidean 
affine space u\(7i) of Hermitian operators with trace 1 - an affine hyperspace of u*(7i). 

It is easy to see that for a compact convex set K in a finite-dimensional affine space E 
the closed convex hull Wn(K°) is just the convex hull con(K°) if only K° C K is closed, 
and that the convex set of positive maps V($(K) is again a compact convex set, this time 
in 21(E). Note that 21(E) is canonically an affine space modeled on the vector space 
2l(E, E) of affine maps from E into E. Moreover, if E is just a vector space, E = E, 
the space %i(E) is a vector space with a canonical decomposition 31(E) = C(E) © E, 
due to the fact that we can write any affine map $ : E — > E uniquely in the form 
$(x) = v(&)(x) + Xq for some v ($) G C(E) and x G E. 

3.1 Fix- extreme positive maps 

In general it is not easy to find extreme points ^(K)° of the convex set of positive 
maps ?P(K), even if extreme points of the convex body K are explicitly known. This is 
exactly the case of the convex bodies ^(T>(1H)) of positive maps in Quantum Mechanics. 

On the other hand, extremality of some positive maps can be established relatively 
easy in the case of maps with many extreme points in the image, as each extreme 
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point in the image fixes partially the map. This is based on the observation that, for 
$ G ty{K), if p G K° is the image p = for some p G K, then <3>j(p) = po for 

any $j of a decomposition $ = ^ Aj$j into a convex combination of G ^(if). 
Indeed, as p$ = = X^^j^Mp) * s a decomposition of the extreme point po into a 
convex combination of points §i(p) G K, then = p$. This immediately implies 

the following. 

Theorem 1. Let K be a compact convex set in an n- dimensional real affine space. If 
a positive map $ G ty(K) has n + 1 affinely independent extreme points in the image 
$(K) of K, then the map $ is extreme positive, $ G ^(K) . 

Proof. Let G X, % — 1, . . . , n + 1, be such that p« = "^(^j) are extreme and affinely 
independent and assume that we have a decomposition $ = t$ + (1 — for certain 
$05^1 G ^P(-?0 and < £ < 1. According to the observation preceding the above 
theorem, $o(?i) = = — P? for alH + 1, . . . , n + 1. But an affine map from 

a n-dimensional affine space is completely determined by its values on n + 1 affinely 
independent points, so $ = $ x = $. □ 

The extreme positive maps $ described in the above theorem (with n + 1 affinely 
independent extreme points in the image ®(K)) will be called fix-extreme positive maps. 

Corollary 1. For any convex body K a positive map ^ which has all extreme points in 
^(K) is extreme positive. In particular, the identity map is always an extreme positive 
map. 

3.2 Example: the closed unit ball in W 1 

Theorem 2. Fix- extreme positive maps of unit balls in Euclidean vector spaces are 
orthogonal transformations. 

The proof of the above theorem will be based on the following lemma. 
Lemma 1. If the function F : M. n — > R, 



where on > and pi G M., i = 1, . . . ,n, has on the unit sphere S*™ -1 = {x G W 
J2i x 1 — 1} local maxima at some n + 1 affinely independent points qi, . . . , q n+ i G S n ~ 
then F is constant on S 12 " 1 . In particular, p\ — ■ • • — p n — and a± — a-i — • • ■ — a. 

Proof. We will use the method of Lagrange multipliers and consider the function 



n 




i=l 
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Since qj, j = 1, . . . , n + 1, are critical points of F, when restricted to the sphere, the 
coordinates (xi, . . . , x n ) of each qj solve the system of equations 

dF\ 

(4) ~dx r ^ = ^ ~~ +Piai = °' i = 1 > ■ ■ ■ ' n ' 

n 



£*, 2 = i- 



i=i 



Moreover, as at qj we have local maxima, the second derivative of F\ must be non- 
positive definite that yields af — A < for all %. We can assume that a± > a<i > ■ ■ ■ > a n , 



so A > a\. 

Assume first that A > a\. Then, 

,-\ _ P%0% 

1 ' Xl ~ A^f ' 

and 



should be equal to 1. But the function G(X) is monotone with respect to of < A < +oo, 
so there is at most one solution (x, A ) of (J4]) with A > o\. It must be therefore at 
least n additional solutions with A = a\. Let k be the number of the biggest «j, i.e. 
ai = a 2 = ■ ■ ■ = Oik > oik+i- We get easily from (jlj) that 

(6) b 1 = --- = b k = 0. 

Hence if k — n then F is constantly a\ on the sphere. It suffices to show now that 
k < n is not possible. Indeed, if x is a solution of (HI) with A = of, we get 

_ PjOtj . _ 

X{ n n i 1 "T~ -Lj • • • J Tb 1 

o{ — ai 

so we can have, together with the solution §5§ corresponding to A > of, an affinely 
independent set of n + 1 solutions only if k — n — 1. But then, due to ([6]), the solution 
((5]) must be 

(7) x = (0,...,0,sgn(p n )), 
so p n ^ and G(A) reduces to 

n( w - iPn a n) 

(A - «2)2 



and 

(8) 
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On the other hand, the solutions corresponding to A = ot\ must be of the form 

VnOtn 



X\ , . . . , X n —\ , „ 

a{ - ai 



so that we have solutions additional to ([7]) only if 

(9) { T n) L < 1 ■ 

But, as Ao > a\i the latter contradicts ([8]): 

(Pn«n) 2 < (Pn« n ) 2 ^ 



(A - a 2 n ) 2 (af - a£ 

□ 



Now we can prove Theorem [2j 



Proof. Let B be the unit ball in an n-dimensional Euclidean vector space E. Let us 
make an identification of E with W 1 with the standard Euclidean norm 

n 

IMI 2 = J2 x2 i- 

i=l 

Let $ be a fix-extreme positive map of B, i.e. $ : IR n — > M. n is an affine map, $(a;) = 
A(x) +p', where A e C(R n ) is a linear map of W 1 such that $(B) C B and $(B) has 
n + 1 affinely independent points • • • , $(^ +1 ) in the unit sphere 5 n ~ 1 . Of course, 

we can assume that q[, . . . , q' n+1 are extreme points of B, so they lie on the sphere as 
well. In particular the map (matrix) A is invertible. 

Now we can apply the singular value decomposition to the matrix A in order to write it 
in the form A = 0\ o T o O2, where Oi, O2 are orthogonal matrices and T is a diagonal 
matrix with positive entries a\, . . . a n > on the diagonal. Since we can write 

= A(x) + b' = 1 o T o 2 (x) + b' = Oi(T o 2 (x) + b) , 

where 0\(b) = b', and since the orthogonal maps preserve B and S 1 ™ -1 , the map &o(x) = 
T(x) + b has the same properties as it is a positive map of B, $o(B) C B and $o(B) 
has n + 1 affinely independent points <3>o((Zi), • • • , &o(<ln+i) on the unit sphere S 1 * 1-1 , 
where g,- = 0^" 1 (gj) are points of the sphere, j = 1, . . . , n + 1. This means that the 
function 

n 

F(x) = \\$> {x)\\ 2 = ^{otiXi+Pif , 

8=1 

reduced to the unit sphere takes at qi, . . . ,q n+ i local maxima. Applying the above 
lemma we conclude that 6 = and F is constant on the sphere, so that $0 — T maps 
the unit sphere into the unit sphere. Hence, T = I and $ = 0\ o 2 is orthogonal. □ 

Remark 1. Theorem [2] can be derived also from the results of [23J. 
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3.3 Example: an extreme map on the plane fixing two extreme 
points 

In the present section we want to present a simple example of an extreme map in two 
dimensions which is a bijection on its two extreme points. To this end let us consider 
the function on the interval [1,1], 

(10) m = fi^v ( i±£V /2 + 2 (i^Y f^f + ^ . 



The function / is concave, hence the subset S of the (x, y) plane bounded by its graph 
and the interval [1, 1] is convex. Let us perform a linear transformation of the (x,y) 
plane, 



11) T : (x,y) h-> 



Under this transformation S is transformed into the set bounded by [1, 1] and the graph 
of 

(12) g(x) = \f(-x), 



i.e. 



1 f 1 + x\ 1 1 — x\ ' ( 1 + X \ ' I 1 — X \ 1 1 — X 



(13) 9 W=^— j {—) {—) + 8 - 

Since /(x) > ^(x) for x G [1, 1], we have T(S) C S*. Moreover T is a bijection on two 
extremal points, (x, y) = (—1,0) and (x,y) = (1,0), of S. Observe also that T is an 
extreme mapping in the sense that for an arbitrary a > 1 there is x G [—1,1] such that 
f(x) — otg{x) < 0, i.e. the linear transformation 

(14) T a :(x,y)»(-x,a^ 
does not map S into S. 

The above described properties of / and g can be established by straightforward calcu- 
lations. Below we illustrate them in Figures 1-3. 
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Figure 2: Second derivative of /. 
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Figure 3: Functions f — g (upper curve) and f(x) — 1.01 • g(x) (lower curve) in the 
vicinity of x = — 1. 



4 Positive maps bijective on pure states - a version 
of Wigner's Theorem 

Before we formulate a version of the Wigner's Theorem jll] let us comment on complex 
antilinear and antiunitary maps in a Hilbert space. A map A : TC — > TC we call antilinear 
if A(ax + (3y) = ax + j3y, where a denotes the complex conjugation of a G C. An 
antilinear map U : TC — >• 7i we call antiunitary if (Ux\Uy) = (y\x) for all x, y G TC. The 
adjoint A* of an antilinear map is an antilinear map defined via the identity 

(Ax\y) = (A^y\x) . 

Any linear (antilinear) map A : TC —>■ TC induces a linear (resp., antilinear) map Ma : 
gl(TC) — > gl(TC) which on one-dimensional maps \x)(y\ takes the form MA{\x)(y\) = 
\Ax)(Ay\. For linear A we can easily represent the map Ma as M^ip) = ApA\ while 
with antilinear maps the situation is a little bit more complicated. 

If an orthonormal basis is chosen then in the Hilbert space TC we can define a complex 
conjugation 

C:TC^TC, C 2 = I, cfca i e?J='22a i e i . 

Instead of Cx we will write simply x. It is clear that (x\y) = (y\x). If A is a complex 
linear map then A = A o C is antilinear and vice versa. Since any continuous complex 
linear (antilinear) map A : TC — ► TC is represented by a (possibly infinite) matrix (a^), 
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where A(ei) = Ylj a l e ji a l so the transposition A \— > A T is well defined: 



A T (e t ) = £ 



j 



and we extend it to the whole TC by complex linearity (antilinearity). For linear A the 
adjoint map can be then written as 

At = C o A T o C , 

so that C o AoC = A T for linear Hermitian A = A^ . If A is antilinear then A = Ao C 
is linear, so 

|Ar)(A/| = \Ax)(Ay\ = Ao \x)(y\ o (2)t . 

But, as easily seen, 

IxVyl =Co|a;)(y|oC, 



so that, for Hermitian p, 

(15) M A (p) = 2p T (2)t. 

The Wigner's Theorem (compare with [13]) can be now formulated as follows. 

Theorem 3. Let ij; : V l (H) -> P 1 (7i) 6e a bisection of pure states in a Hilbert space 
TC preserving the transition probabilities 

(16) (V(pi)l^(P2))* = (pi|p 2 )*. 
T/ien, i/iere zs a unitary U : TC — > 7i st/c/i £/ia£ 

(17) ip(p) = U pU^ for all pure states p , 
or 

(18) ^(p) = Up T U^ for all pure states p, 

where p i— > p T zs i/ie transposition associated with a choice of an orthonormal basis in 

n. 

The standard versions of Wigner's Theorem usually consider (unit) vectors of the 
Hilbert space rather than pure states. But if x, y are unit vectors representing pure 
states pi and p 2 , respectively, then 



(pi,P2>* = \(x,y)\ 



2 

> 



so that preserving \ (x, y)\ is the same as preserving (pi,p2)*. Moreover, any unitary (or 
antiunitary) action in the Hilbert space, x \— > Ux, induces on pure states p = \x)(x\ the 
action (fTTl) or (fl8|) . We will call the maps (iTTj) and (fl8l) defined on pure states or on 
u*(TC) Wigner maps. The Wigner maps on u*(TC) can be abstractly characterized as 
follows. 
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Theorem 4. A linear map ip : u*(7i) — > u*(7i) is a Wigner map if and only if it is 
positive and orthogonal. 

Proof. The Wigner maps are clearly positive and orthogonal, so let us assume that if) 
has these properties. For all Hermitian pi,p2 we have therefore ( fl6l) and we know that 
ip(p) is positive semi-definite if p is. The map ip is orthogonal, therefore invertible, 
and its inverse ip~ l is orthogonal as well. Let us observe that ip~ x is also a positive 
map. Take a pure state p and suppose that ip~ l {p) has the spectral decomposition 
vb~ 1 {p) = p + — p_ into a difference of positive semi-definite operators p+,p~ which are 
orthogonal, (p + |p_)* = 0. Then p is a difference of orthogonal positive semi-definite 
operators p = ip(p+) — ip(p-) and, as p is a pure state, ip(p-) (thus p_) must be 
0. A similar argument shows that the image ip(p) of any pure state p is a positive 
semi-definite operator which is not decomposable into a sum of orthogonal positive 
semi-definite operators, so ip(p) is a pure state up to a constant factor. Since 

Tr(^(p) 2 ) = ^(p)|^(p)) = (p|p) = l, 

this factor equals 1 and we conclude that ip induces a bijection on pure states. □ 

We will now prove a theorem which extends Wigner's Theorem and which relates it to 
the problem of extreme positive maps. 

Let u*f(7{) be the linear subspace of u*{7i) consisting of Hermitian finite-rank operators. 
For Ki,K 2 C u\(7i) we say that a map ip : K\ — > K 2 is affine if if) is the restriction to 
K\ of a trace-preserving linear map $ :< K\ >^< K 2 > from the linear span < K\ > 
of K 1 in u*{H) into the linear span < K 2 > of K 2 in u*(Tt), $(-^1) C K 2 . 

Theorem 5. Let ip : T' 1 (7i) — > P 1 (7i) &e a bijective map. The following are equivalent: 

(a) ip is affine; 

(b) ip preserves transition probabilities between pure states; 

(c) ip is a Wigner map. 

If any (or all) of these cases is satisfied, there is a unique continuous affine extension 
* : u\(H) -> u\(H) ofvb which is extreme positive, $ G ^(V(H)f . 

Proof, (a) =^> {b). Since u*JTL) is spanned by the set T) X {J-C) of pure states, let $ : 
u*f(TC) — > u*j-(H) be the (unique) linear trace-preserving map on the space u%(Tl) of 
finite-rank Hermitian operators inducing t/> on T' 1 (7i). Since $ maps convex combi- 
nations into convex combinations, $ maps finite-rank density states into finite-rank 
density states, so $ is a positive map. We will prove that $ is a linear isomorphism on 
u*f(l~C). This follows from the fact that $ preserves the rank, i.e. induces a bijection on 
V k {H) for each k = 1,2,.... 

To see the latter let us remark that the rank, rank(p), of p G u*j{7i) is defined as a 
minimal number of pure states whose linear combination is p, so that, as $ is linear 
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and is a bijection on pure states, rank(&(p)) < rank(p). Conversely, if p G T>{Ti) 
is of rank k, then it is a convex combination of some pure states pi, . . . ,pt, thus the 
image by $ of a convex combination of pure states $ _1 (pi), . . . , $ _1 (pfc). This shows 
that T> k (7i) C ^{T> k {7i)). The linear map $ induces a bijection on 'D 1 (7Y), so assume 
inductively that it induces a bijection on T> l (Ti) for / < /c. Let now p G u*f(Ti) be of 
rank &; + 1, p G V k+1 (H), with the spectral decomposition p = Xli=o ^P«> with A« > 0, 
^ Aj = 1, and p« = |a;j)(xj| being pairwise orthogonal pure states, (xi\xj) = 0, % ^ j. 
We must show that the rank of <J>(p) is k + 1. 

Suppose the contrary. Hence, according to the inductive assumption, p' = $(p) is of 
rank k. As p = A po + P, where p = X^=i A^pj, thus p 7 = $(p) is of rank k, the image of 
$(p ) = l^oX^ol? thus x' belongs to the image of p 7 . Consider the spectral decomposition 
p 7 = X^i=i Kp'ii with A 7 > 0, J2 i A 7 = 1, and p\ = |^)(^| being pairwise orthogonal pure 
states, (a^Xj) = 0, i ^ j. Let 7i (resp., 7^q) be the subspace in 7i spanned by the 
vectors xi, . . . ,x k (resp., x' 1: . . . ,x' k ) and let V l (Ho) (resp., V 1 (H' )) be the set of all 
pure states of Ho (resp., Hq), i.e. these pure states from V l {7i) which are represented 
by unit vectors from Tio (resp., 7~C' ). It is now clear that p' = $(po) G 'D 1 (7io). Note 
that pure states rj from 'D 1 (7io) can be characterized as such pure states which added 
to p 7 do not change the rank, rank(p r + rj) — rank(p') = k. According to the inductive 
assumption this implies that rank(p + $ _1 (?7)) = k as well, but 

rank(p + $ _1 (p 7 )) = rank(p + p ) = rank(p) — k + 1 , 

which is a contradiction. 

Since we know now that <£> induces bijections on each V h (H), k — 1, 2, . . . , it is easy 
to conclude that it is a rank-preserving isomorphism, so that is also a positive 
map. Indeed, as P 1 (7i) spans u*j(7i), it is clearly "onto". It is also injective, since 
$(Ap — A 7 p 7 ) = 0, where A, A 7 > and p, p 7 G T>(1~L) are density states of finite ranks, 
implies ($ is positive) that $(p) = $(p 7 ) = 0, thus p = p 7 = 0, as $ preserves the rank 
of density states. 

To finish the proof, we will need the following lemma. 

Lemma 2. Let p G V(TL) be a density state of a finite rank k. Then the square of the 
Hilbert- Schmidt norm \\p\\l = Tr(p 2 ) can be characterized as the maximum of the expres- 
sions Yli=i over a M decompositions p = Yli=i ^iPi of p as a convex combination of k 
pure states pi,p2, • • • ,Pk- This maximum is associated with the spectral decomposition. 

Proof. As 

(19) Tr(p 2 ) = J2 X " + 2 J2 M/toM' 

i i^j 

and 



(20) 



A i A j Tr(p i p i ) > , 
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we have 

k 

(2i) ll/€>£ A ? 



i=l 



Moreover, we have equality in ( 12TI) if and only if we have equality in ( l20i) . so as all 
Aj > 0, if and only if (pi\pj)* = Tr(p;Pj) = for all i ^ j. This corresponds to the 
spectral decomposition. □ 

The above lemma implies that the map $ preserves the Hilbert-Schmidt norm of density 
states. Indeed, if we use the spectral decomposition to write p as a convex combination 
p = $^»=i AjPj of pure states, then $(p) can be expressed as a convex combination 
of pure states $(p) = Yli=i ^i^(Pi) w hh the same coefficients, so that ||$(p)||* > 
Ylii ^1 = IIpII*- we can apply the above consideration to $ _1 instead of $ and get 
ll^ -1 ^) II* — IMI* f° r an y density state of rank k, in particular for 77 = $(p). We get 
therefore \\$(p)\\l = \\p\\l 

Let us now take two pure states pi, P2 and consider p = \{p\ + P2)- Since, according to 



and 



we have 



thus 



Ml* = o i 1 + (Pi\pi)< 



\p\\1 = \Mp)\\1 = \\1(Hpi) + Hpi))\\1, 



\{1 + {P1\P2).) = \(1 + (*(P1)\*(P2))*) 



(pi|p 2 >* = ($(pi)|$(p 2 )>* = (1P{PI)\1P{P2)}*, 

so ip preserves transition probabilities between pure states. 

(b) => (c) is the Wigner's Theorem. 

(c) =>• (a) is obvious. 

Moreover, if) has an obvious unique continuous extension \1/ : u\(l-L) — > u^(7i), ^(p) = 
[7 pC/' or ^(p) = Up T U^ . Since \l/ is positive and has all extreme points in its image, 
it is extreme positive according to the obvious infinite-dimensional version of Corollary 

m □ 

If the dimension of the Hilbert space Ti, is n, then the dimension of the affine space 
m*(7Y) of Hermitian operators with trace 1 equals n 2 — 1 and we know from the general 
theory (Theorem CD) that a positive map $ G tyiViTi)) possessing n 2 pure states in the 
image $(22(7^)) of T>{7i) is extreme positive (we called such positive maps fix-extreme). 
We finish this section with the following conjecture motivated by Theorem [2j 

Conjecture 1. Any fix- extreme positive map $ G ^P(P(7^))° is a Wigner map. 
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5 Completely positive maps acting on pure states 

A linear map A : u*(H) — ► u*(H) is called completely positive (CP) if A® I N : u*(H) ® 
M^r ~~ ► u*(7i) ® Mjv, where M^r is the algebra of complex N x N matrices, is positivity- 
preserving for all N. It was shown by Choi [21] and Kraus [25] that each CP map 
admits a representation in the so called Kraus form 

s 

(22) Ap = ]T Vtfpis, 

i=l 

where V k are operators acting on 7i. 

Definition 2. We say that a CP operator ( 1221) acting on Hermitian operators of a 
Hilbert space 7^ is extreme if any decomposition, A = A\ + • • • + A r , into a sum of CP 
operators A\, . . . , A r is irrelevant, i.e. all the operators Ai, . . . , A r are proportional to 
A, Ak = a^A, for some G R + . In particular, all operators Vi are proportional and A 
can be written as a single Kraus operator 

Ap = V ] pV. 

Not to deal with the cone of CP operators but rather with a compact convex set (if the 
dimension of 7i is finite), one has to put certain normalization condition. We can try 
to use, for instance, the trace Tr(A) of a CP operator (1221) as a linear operator on the 
complex vector space gl(7i) which is the same as the trace of (1221) as an operator on 
the real vector space u*(7~C) of Hermitian operators. 

Theorem 6. If A is a CP operator in the form then 

(23) Tr(A) = ^|Tr(^)| 2 . 

i=i 

Proof. It suffices to prove (1231) for a single Kraus map A = My- Choose an orthonormal 
basis {e{\ in Ti and write V in this basis as a complex matrix V = (%•). As an 
orthonormal basis in gl{7i) we can take pjk = |ej)(efc|. We have 

Tr(M v ) = ^{p jk \V^p jh V) = Y t Tr{p kj V^p jh V) 

= ^(eklpkjV^pjkVek) = ^ VjjV kk = |Tr(V)| 2 . 

j,k j,k 

□ 

As we can see, the trace can vanish for non-zero CP operators which makes the nor- 
malization by trace impossible. There is, however, another possibility of normalizing 
CP operators provided by the Jamiolkowski isomorphism [26] 



(24) 



J : gl(gl(H)) - gl(gl(H)) 
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The Jamiolkowski isomorphism maps the operator M^(p) = ApB^ into the rank-one 
operator In particular, CP operators correspond, via the Jamiolkowski isomor- 

phism J, to positive semi-definite operators on gl(TC) (see, for instance, |20j). Any 
Kraus operator My(p) = VpV^ corresponds to the one dimensional Hermitian operator 
|V)(V|. The spectral decomposition of J {A) results in the decomposition (|22|) with 
Vi being mutually orthogonal with respect to the Hilbert-Schmidt product Tr(A'B) in 
gl(TC). Such a decomposition of a CP operator we will call a spectral decomposition. We 
will call a CP operator normalized if it corresponds, via the Jamiolkowski isomorphism, 
to a trace-one operator. If fl22|) is a spectral decomposition of a CP operator, then it is 
normalized if 



We will call Tr s (v4) = Tr ( Yli=x ^« ) the spectral trace of a CP operator A. The 
spectral trace is greater or equal to and equals only if A = 0. The convex set of 
normalized CP operators on gl{7i) we will denote by NCP(7i). It is clear that extreme 
points in NCP{H) are exactly single normalized Kraus maps. 

From now on we assume the Hilbert space Ti, to be of a finite dimension n. 

Recall that on the space of Hermitian operators on 7i we have a canonical scalar product 
(p, p% = Trpp' and that V(H) denotes the cone of positive semi-definite operators. We 
will call elements p\, . . . , pk, k > n, of V{7i) to be in general position if for any non-zero 
p G V(H) we cannot find n of them which are orthogonal to p. If pi are of rank-one, 
Pi = l^iX^il, where | Xj) G H, this simply means that any n vectors of | x%), ■ ■ ■ , \ Xk) 
form a linear basis in Ti. 

Theorem 7. The following statements are equivalent: 

(a) A is invertible and A^ 1 is a CP operator. 

(b) A is invertible and extreme, i.e. 



with V - invertible. 

(c) For any set of pure states p±, . . . , p n+ i in general position, Ap%, . . . , Ap n+ i are 
operators of rank one in general position. 

(d) In the image A(P(7i)) there are n + 1 rank-one operators in general position. 

(e) There exists a set of pure states pi, ... , p n+ i in general position such that Ap%, . . . , Ap n+ i 

are operators of rank one in general position. 

Proof. (a)=^(b). It was proven in [27] (Theorem 7). 

(b)=^(c). If p = \x)(x\ then V^pV = \V^x)(V'x\, hence A maps pure states into rank-one 
operators, and V, being invertible, preserves all linear independencies. 




Ap = V j pV 
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(c) =>(d) - trivial. 

(d) =r-(e). Let pi, ■ ■ ■ , p n +i G P(7~0 be positive semi-definite operators such that rjj = 
A-Pji j = 1, • • • ,n + 1, are rank-one operators in general position. First, let us remark 
that we can assume that pj are density states, since proportionality plays no role here. 
Second, we can assume further that they are pure. Indeed, if p is a state with the 
spectral decomposition p = J^ fc a k ^ k into a convex combination of pure states and 
if Ap = i] is rank-one positive semi-definite operator, then r\ = Ap = J2 k a k A^ k is a 
convex combination of rank-one positive semi-definite operators A£ k , so r\ is positively 
proportional to A£ k for all k, since pure states are extreme points in the convex body 
of all states. By a similar argument, all states V} pj V$, i = 1, . . . , s, are proportional to 
r)j, say V- pjVi = (3jT)j, where, of course, YH=iPi = 1- 

Let us write pj = \xj)(xj\, rjj = \yj)(yj\, j — 1, . . . , n + 1, for some vectors | x-,-), | y^) . We 
claim that the states pi, . . . , p n +i are in general position, i.e. any n among the vectors 
| Xi), . . . , | x n+ i) are linearly independent. 

For, assume the contrary, i.e. that, say | xi), . . . , | x n ), are linearly dependent. Hence, 
one of the vectors, say | x n ), can be written as a linear combination | x n ) = b±\ x±) + 
• • -b n -i\ x n -i). Since V} pjVi = \ V^Xj){V^Xj\ is proportional to rjj = \yj){yj\, the vector 
| y}xj) is proportional to \yj) for j = 1, . . . , n + 1. In particular, all the vectors | V}xj) 
with i — 1, . . . , s and j — 1, . . . , n — 1, belong to the linear span span(| yi), . . . , | y n -i)) 
of the vectors | y± ),..., | y n -i)- Therefore 

^/l ^n) = fell? | ^i) H h K-iVi\ x n _i) G span(| yi), . . . , | y„_i)). 

Since ^|a; n ) is proportional to |y n ), and | yi), . . . , | y n ) are linearly independent, we 
have Vi | x n ) = for alH = 1, . . . , s. But this contradicts the property Ap n = r] n , i.e. 

T,t=l V i\ X n) = \Vn)- 

(e) =^(a). Put pj = \xj)(xj\, rjj = Apj = \yj)(yj\ for some | Xj), \ yj) e Ti, j — 1, . . . ,n+l. 
Since rjj = Apj is of rank one, all the positive semi-definite operators V} PjV% of its 
decomposition 

s 

i=i 

must be proportional to rjj, V} PjVi ~ rjj. This, in turn, means that V}\ Xj) are propor- 
tional to | yj), 

V i\ x o) = a i\yj), i = l,-..,s, j = l,...,n + l. 

As any n vectors of | x±), . . . , | x n +i) are linearly independent, all the coefficients dj of 
the decomposition 

| x n+1 ) = ai| xi) H h a n \ x n ) 

are non-zero, a.j ^ 0. Since 

\//| rc n+ i) = aiV/| xi) H h a„\//| rr n ) = aia]\ y\) H h a„a™| y n ) 

are proportional to | y„+i) and | y n +i) has a decomposition | y n+ i) — b\ \ yi) + - • - + b n \ y n ) 
into a linear combination of linearly independent | yi ),..., | y„), the vector (a], . . . , a") G 
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C n must be proportional to (aj/fei, . . . , a n /b n ) G C n . In consequence, it means that the 
operators V- are proportional to the operator uniquely defined by the conditions 

V^\x J )= b My 3 ). 

Hence 

Ap = (3V^pV 

for some (3 G 1R+ and, clearly, (3 ^ 0. □ 

Remark 2. It is not enough to take only n pure states pi, ■ ■ ■ , p n - Let us take, for 
example, Vi diagonal, Vi = diag (Xj , . . . , A"), ^ 0, and not proportional. If now 
Xj = Xf, we have an infinite set of pure states that are mapped to rank-one operators 
and any n of them are not in general position, but A = Y2i=i V} P^i ^ s n °t extreme. 

Corollary 2. If a CP operator 

s 

(25) Ap = ^ V?pVi 

i=i 

is trace-preserving (resp., unity-preserving) and such that the image of density states 
A(T>) contains n + 1 pure states in general position, then A is unitary, Ap = WpU. 

Proof. It is enough to make use of Theorem [7] and observe that trace-preserving (unity- 
preserving) yields VV^ = I (resp., V*V — I), so V is unitary. 

□ 



6 Extreme bistochastic CP maps 

Extreme unity-preserving CP maps have been described by Man-Duen Choi [23]. We 
can reformulate his result as follows. 

Theorem 8. Let NCPi{ri) (resp., NCP^iTi) ) be the convex body of normalized unity- 
preserving (resp. trace-preserving) CP maps on gl{H). Then A G NCPj(7i) (resp. 
A G NCP^Ti)), with the spectral decomposition 

s 

(26) Ap = ^ vfpVt 

i=l 

is extreme if and only if the operators {V^Vj : i,j = 1, . . . , s} (resp., {ViVj : i,j = 
1, . . . , s}) are linearly independent in gl(TC). 

Proof. We will sketch a proof making use of the Jamiolkowski isomorphism f l24|) which 
associates with the CP map fl26|) the Hermitian operator 

s 
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on gl(7i). Extreme normalized CP operators correspond therefore, via the Jamiolkowski 
isomorphism, to pure states on the Hilbert space gl(TC). On the space u*(gl(TC)) of Her- 
mitian operators acting on the Hilbert space gl(H), in turn, we can define two canoni- 
cal M-linear maps F^,F 2 : u*(gl(Tl)) — > u*(7i) which associate with rank-one operators 
\V)(V\ the operators F^V^V |) = and F 2 {\V)(V\) = VW, respectively. The unity- 
preserving CP operators A correspond therefore to Hermitian operators constrained by 
the equations F\(J(A)) = I. The corresponding extreme points J(A) e u*(gl(H)) 
need not be pure states. It suffices that the level sets of / of the linear constraint F\ 
are transversal to the face of the point, i.e. the function F\ has the trivial kernel on 
the tangent space Tj(A) of the face at J {A). But this tangent space is known (see 
e.g. [3 [27]) to consists of all Hermitian operators with the range equal to the range of 
J{A\ i.e. the operators of the form Ylij ^l^i X^j Ij "where (A u ) is a Hermitian matrix. 
Hence A is extremal in NCPj(7i) if and only if 

(27) iME^Wx^i) = E A ^ f ^ * G < n ) 

for all Hermitian (A y ) ^ 0. As we can decompose any operator into the sum of a 
Hermitian and an antihermitian ones, we can rewrite (1271) in gl{7i) instead of u*{7i) 
using an arbitrary complex matrix (A w ) 7^ 0. This is nothing but the complex linear 
independence of the operators {V^Vj : i,j = 1, . . . , s} in gl(H). For trace-preserving 
CP operators the reasoning is identical with the function F 2 replacing F±. □ 

The maps from NCP*(H) = NCP^H) H NCP Tr (H) are sometimes called bistochas- 
tic. The above understanding of the Choi's result gives us easily a characterization of 
extreme bistochastic maps. 

Theorem 9. A bistochastic map A e NCP*(1~C) with the spectral decomposition 

s 

Ap = J2 V}pVt 

i=l 

is extreme if and only if the operators {V^Vj © ViVj : i,j = l,...,s} are linearly 
independent in gl(H © H). 

Proof. The proof is analogous to the above one with the difference that our constraint 
function is now 

(F 1; F 2 ) : u*{gl{H)) -> u*{H) x u*{H) ~ u*{H) © u*{H) . 

The condition (I27]) is therefore replaced by the condition 

\ ij V}Vj ^ or X'ViV} ^ for all Hermitian (A iJ ) ^ . 

We can rewrite it as 

X>"(\y\) ) / 

and pass to an arbitrary complex (A 4 - 7 ) 7^ like before. □ 
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7 Examples 



Let us illustrate some of the previous reasonings and results in the simplest cases of 
maps on states on two- and three-dimensional spaces. In the following subsection we 
show three examples of extreme maps: a generic one possessing exactly two pure states 
in its image, a non-generic one with only one pure state in the image, and an extreme 
map having a continuous family of pure states in the image. The last one, according to 
Theorem [2j is not completely positive but a merely positive extreme map. Note that 
these cases have been considered also in [291. 



In the second subsection we give an example of an extreme completely positive map 
acting on C 3x3 which does not have any pure state in its image. Such a situation is 
impossible for maps acting on qubits (i.e. maps on C 2x2 ). 

7.1 Extreme completely positive, positive, stochastic and bis- 
tochastic maps for n = 2 

A state on C 2 can be parameterized by a unit vector (x, y, z) e M 3 



(28) 

where 

(29) 



p = - (I + xg x + yo 2 + za 3 ) 



0"! 



" 


1 " 




' -i ' 




" 1 


1 





, o 2 = 


i 


, ^3 = 


-1 



A positive trace-preserving map 

(30) p ^ p' = - (I + arVj + y'a 2 + z'a 3 ) 

is thus determined (up to rotations which are irrelevant from the point of view of this 
paper) by four parameters Ai, A 2 , A 3 , t, such that 



(31) 



x' = Xix, y' = X 2 y, z' = X 3 z + 1. 



The parameters Ai, A 2 , A 3 ,t must fulfil particular conditions to ensure the positivity of 
the map (see [28], [29] for details). 

The image of the unit sphere x 2 + y 2 + z 2 = 1 under (130]) is the ellipsoid 



(32) 



+ 



z-t 



i. 



Obviously, for a positive map (13"U1) the ellipsoid is inside the unit sphere. For extreme 
maps it has to have points on its surface common with the surface of the unit sphere (i.e. 
some pure states are mapped into pure states). For extreme CP maps two possibilities 
occur [291: 
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1. Ai = cos(w), A2 = cos(f), A3 = cos(u) cos(v), t = sin(w) sin(t>), < < u < v. In 
this case the ellipsoid (1521) has three different axes and it touches the unit sphere 
at two points (see Fig. @|. 




Figure 4: An extreme CP map having exactly two pure states in the image. Projections 
along two perpendicular axes. 



2. Ai = A2 = cos(w), A3 = cos 2 (w), t = sin 2 (w), in which case the ellipsoid ( |32l) 
touches the unit sphere at a single point x = y = 0, z = 1. (see Fig. [5]). 




Figure 5: An extreme CP map having exactly one pure state in the image. Projections 
along two perpendicular axes. 



3. Geometrically it is obvious that without upsetting the extremality of the map we 
can make the ellipsoid ( 1321) touching the unit sphere along a full circle (see Fig. [6]). 
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In this case 

Ai = A 2 = a/1 — cos 2 (m) cos 2 (t>) , 
A3 = sin (it) 1 — cos 2 (w) cos 2 (t>) , 
t = sm(u) sin 2 (t>) . 

For n^O^w the map is definitely not a unitary one (its image is a proper subset 
of the unit sphere) and in its image there are more than 3 pure states (in fact 
the whole circle of states at which the ellipsoid touches the unit sphere). From 
Theorem [2] it follows thus that the map cannot be a completely positive one. 
Indeed, for the chosen values of Ai, A2, A3, and t the map is an extreme positive 
[23] . The fact that it is not completely positive can be checked independently by 
finding that its image under the Jamiolkowski isomorphism [26] is not positive 
semi-definite. 




Figure 6: An extreme positive map having a continuous family of pure states in its 
image. Projections along two perpendicular axes. 



7.2 An extreme completely positive map having no pure states 
in its image 

Let us consider a CP map on C 3x3 defined by the following Kraus operators 



(33) 
Vi = 



l/y/3 
l/y/2 















L VT+~c 



l/y/3 

l/y/2 





l/y/3 






A straightforward calculation gives ViV-/ + V2V2 + V3V3 = I. 
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For a = the matrices ViYj form a basis in the space of 3 x 3 matrices, hence it is also 
true for small a. The map Ap = Yli=i ^1 pVi ls thus an extreme CP map (theorem EJ) . 
For q ^ there is no | y) such that V±\ x) ~ | y) for some | x) and i — 1,2, 3, hence A 
does not send any pure state into a pure one. 
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